Gamma-ray background anisotropy from galactic dark matter substructure 
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Dark matter annihilation in galactic substructure would imprint characteristic angular signatures 
on the all-sky map of the diffuse gamma-ray background. We study the gamma-ray background 
anisotropy due to the subhalos and discuss detectability at Fermi Gamma-ray Space Telescope. In 
contrast to earlier work that relies on simulated all-sky maps, we derive analytic formulae that 
enable to directly compute the angular power spectrum, given parameters of subhalos such as mass 
function and radial profile of gamma-ray luminosity. As our fiducial subhalo models, we adopt 
M~^'^ mass spectrum, subhalos radial distribution suppressed towards the galactic center, and 
luminosity profile of each subhalo dominated by its smooth component. We find that, for multipole 
regime corresponding to 6' < 5°, the angular power spectrum is dominated by a noise-like term, with 
suppression due to internal structure of relevant subhalos. If the mass spectrum extends down to 
Earth-mass scale, then the subhalos would be detected in the anisotropy with Fermi at angular scales 
of ~10°, if their contribution to the gamma-ray background is larger than ~20%. If the minimum 
mass is around IO'^Mq, on the other hand, the relevant angular scale for detection is ^1°, and the 
anisotropy detection requires that the subhalo contribution to the gamma-ray background intensity 
is only ~4%. These can be achieved with a modest boost for particle physics parameters. We also 
find that the anisotropy analysis could be a more sensitive probe for the subhalos than individual 
detection. We also study dependence on model parameters, where we reach the similar conclusions 
for all the models investigated. The analytic approach should be very useful when Fermi data are 
analyzed and the obtained angular power spectrum is interpreted in terms of subhalo models. 

PACS numbers: 95.35.-l-d, 95.85.Pw, 98.35.Gi, 98.70.Vc 



I. INTRODUCTION 

Modern astrophysical and cosmological measurements 
strongly support existence of nonbaryonic dark matter. 
Although the true identity of dark matter is unknown ob- 
servationally and experimentally, there are several well- 
motivated particle-physics models that provide a candi- 
date particle for dark matter. Weakly interacting massive 
particles (WIMPs) such as supersymmetric neutralinos 
are perhaps the most popular candidate P, [1] . As they 
interact with themselves as well as standard-model parti- 
cles, many experiments are being carried out to look for 
signatures of scattering of WIMPs off nuclei in under- 
ground detectors and of self-annihilation of WIMP par- 
ticles in dark matter halos 

Recently launched Fermi Gamma-ray Space Tele- 
scope Q features promising capability to detect gamma 
rays from WIMP annihilation in the right energy 
range ,5]. In addition, recent numerical simulations 
find that dark matter in a host halo is distributed 
in clumpy substructure (subhalos) [1, S B 111 whose 
masses raiige quite widely, potentially down to Earth 
mass im, in, in, [l3|- This feature is encourag- 
ing because annihilation probability is proportional to 
density squared, and thus the significant dumpiness 
boosts the gamma-ray yields. Several avenues have 
been proposed for Fermi to look for annihilation gamma 
rays: galactic center [l^, EB; IHj E^, relatively 
massive substructure often associated with dwarf galax- 
ies 20, 21, 22, .23, ^5, 26, 27, 28, 29, 30, 31], proper 
motions of nearby small subhal os (3p. l33l . d iffuse gamma- 
ray background 0, [H, [H [13, ialaEoi , etc. 



There is an increasing interest in statistical analysis 
of the all-sky map of the gamma-ray background ob- 
tained with Fermi in the near future. References 0,111 
computed angular power spectrum of the gamma-ray 
background from annihilation in extragalactic dark mat- 
ter halos, and showed that the signature would be dif- 
ferent from that of ordinary astrophysical sources (see 
also, Refs. [H, HI, [H, H^]). The same approach has 
been applied to signals from substructure in the galac- 
tic halo 47,,^]. This is indeed important, because the 
galactic substructure typically gives a larger contribution 
to the diffuse gamma-ray background than extragalactic 
halos do m, 111, 113 . In Refs. [13, 13, first the all-sky 
gamma-ray map was simulated from sets of subhalo mod- 
els, and then the map was analyzed to obtain the angular 
power spectrum. In addition, one-point probability dis- 
tribution function of the gamma-ray flux has also been 
studied 

In this paper, we revisit the gamma-ray background 
anisotropy from dark matter annihilation in the galactic 
subhalos. In contrast to the earlier works [13, 13 that 
heavily relied on mock gamma-ray maps generated from 
subhalo models, we develop an analytic approach to com- 
pute the angular power spectrum directly. This way, we 
are able to calculate the angular power spectrum easily 
and more quickly, if we specify some input parameters 
and characteristics of galactic subhalos. This would be 
in particular useful when we have results of actual Fermi 
data analysis, and try to give physical interpretation for 
them. 

We find that the angular power spectrum Cg is divided 
into two parts: one depending on (ensemble-averaged) 
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distribution of subhalos in a host ( "two-subhalo" term, 
Cg^^^), and the other depending on the emissivity profile 
of single subhalos as well as the number of subhalos sig- 
nificantly contributing to the background intensity ( "onc- 
subhalo" term, Cj^^^). The latter would be shot noise if 
the subhalos were completely point sources, but Fermi 
will be able to see deviations from the shot noise due 
to the angular extension of the relevant subhalos. Using 
the latest subhalo models following recent numerical sim- 
ulations, we give predictions for the angular power spec- 
trum, and show that the multipole range 10 < £ < 100 
would be a favorable window for anisotropy detection. 
We also discuss the detectability of the angular power 
spectrum from subhalos with Fermi, which turns out to 
be promising, potentially better than the detection of 
subhalos as identified gamma-ray sources. 

This paper is organized as follows. In Sec. |TT1 we give 
formulation for angular power spectrum as well as mean 
intensity of the gamma-ray background. These formu- 
lae derived are applied to several subhalo models in the 
subsequent sections. In Sec. IIIIl we study a simple case 
in which all the subhalos are assumed to be a point-like 
gamma-ray emitters. The case of extended subhalos is 
addressed in Sec. llVi where we also discuss detectability 
with Fermi. We close this paper by discussing the results 
in Sec. IVland by giving concluding remarks in Sec. IVII 



II. FORMULATION 

A. Relevant quantities of subhalos 

We assume that a fraction / of the mass of the galactic 
halo is in the form of subhalos, and 1 — / is distributed 
as a smooth halo. For the density profile of the Milky- 
Way dark matter halo, we adopt a spherically symmetric 
Navarro-Frenk- White (NFW) profile [53|: 



Pmw,nfw('") 



Ps,MW 



(r/rs,Mw)(l + r/rsMw) 



(1) 



where r is the galactocentric radius, r^^MW and Ps,MW 
are the scale radius and scale density of the Milky- Way 
halo, respectively. This profile extends up to a virial 
radius rvir,MWj and an enclosed mass within this radius 
is defined as a virial mass Mvir.MW- We use the following 
values for these parameters: r-g^MW — 21.5 kpc, ps,MW = 
4.9 X lO^Mfp kpc~^, rvir mw = 258 kpc, and Mvir mw — 
IQi^M© Hf. 

As subhalo number density per unit mass range, we 
define a subhalo mass function, dushir, M)/dM, and 
upper and lower limits of the function by Mmax and 
-^minj respectively. Numerical simulations imply that 
the shape of mass distribution follows typically a power 
law, which is close to dush/dM oc M~^, i.e._jthe same 
amount of subhalo masses per decade [H, [H, [13, [H, [s^ . 
We also assume that there is a one-to-one relation be- 
tween subhalo masses and luminosities L(M), and there- 
fore, the luminosity function is written as dush/dL = 



{dush/ dM)\dM / dL\. The upper and lower limits on the 
luminosity function are then given by Lmax = ^(.^max) 
and Lnun = i(-^'^min), respectively. After integrating the 
mass (luminosity) function over mass (luminosity), we 
obtain the number density of subhalos rtgh (f) ■ 

We assume that each subhalo has extended, isotropic 
emissivity profile around its center (we call this "seed" 
position), u(rsh, M)L, with the profile function u(rsh, M) 
normalized so that it gives unity after volume integra- 
tion. We also define the Fourier transform of u{rsh,M): 
u{k, M), where k is the wave number. 

The most relevant equations are Eqs. ([5]) and (|18 p ~ ([20|) 
derived in the remainder of this section. The readers who 
are only interested in application of these equations to 
subhalo models may skip to Sec. IIIIl 



B. Gamma-ray intensity from subhalos 

We label positions of seed of a subhalo i hy Xi, and its 
luminosity and mass by Li and Mi, respectively. With 
these definitions, the gamma-ray intensity towards a di- 
rection n is given by the line-of-sight integration (ds) of 
the emissivity: 



/(n) = 



^ J d.s'Y^uisn - Xi,Mi)Li 

i 

= ^ J ds J dL J d^x^5^{x - Xi)5{L - Li) 



X u{sii — X, M)L, 



(2) 



for one realization of the universe, where 6^ is the A^- 
dimensional delta function. Throughout this paper, we 
define the intensity as a number of gamma-ray photons 
per unit area, time, and solid angle, and the luminosity 
as a number of photons emitted per unit time. We also 
assume > 10 GeV as a targeted gamma-ray energy. 

We now take ensemble average over infinite number of 
realizations of the universe. The discrete source distribu- 
tion then becomes continuous function; i.e., 

l^5\x^x.)5{L-U)'^^''^^^^, (3) 

where the bracket represents the ensemble average. Using 
these in Eq. an ensemble-averaged intensity is 



X u(sn — X, M)L. 



(4) 



We further assume that spatial extension of each sub- 
halo is much smaller than a scale on which subhalo dis- 
tribution significantly changes. With this reasonable as- 
sumption, we could take the luminosity function dugh/dL 
out of the volume integral by taking x w sn, since it is 
a slowly varying function of a;. As the integration of 
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u{sn — X, M) over x simply becomes unity, the ensemble 
average of intensity is 



1 

47r 



dL 



.{L) 



dL 



. (5) 

where we specified upper and lower limits of the integrals. 
Galactocentric radius corresponding to sn (appearing as 
index of luminosity function) is obtained through the re- 
lation: = Tq + s'^ — 2rQS cosip, where Tq — 8.5 kpc 
is the galactocentric radius of the solar system and ip 
is the angle between n and the direction to the galac- 
tic center. We obtain the lower limit of the s-integral 
by the detection criterion L — AnslFscns with the flux 
sensitivity of Fermi (typically Fgcns — 2 x 10^^° cm~^ 
s~^ for photons that we consider). By setting this, we 
do not add contributions from subhalos bright enough 
to be identified as individual sources. The upper limit 
Smax(iT-) corresponds to Tvir.MW through the relation 



' vir,MW 



-5 max cosV'- 



By further averaging over directions n, we obtain a 
mean gamma-ray intensity 



1 



a 



sky 



1 

47r 



— / dL 



ds^^L, (6) 



where we represent quantities averaged over n by putting 
a horizontal line on top of them, and ilsky is the solid 
angle of the sky over which the averages are taken. In 
the second equality, we approximate Smax ~ ''vir,MW, as 
7*0 ^ 7'vir,Mw; since Smax is uow independent of n, we 
could perform the angle average of dugh/dL before inte- 
grating over s and L. 

The expected number of subhalos that would be de- 
tected as individual sources from figky is then 



(^sh) — ^^sky 



dL 



s,{L) 



dss 



tnshjs, L) 
dL 



(7) 



as we count sources close enough to give a flux larger 
than Fscns; i.e., s < s*(L). 



C. Angular power spectrum of gamma-ray 
background from subhalos 

We decompose the gamma-ray intensity map with 
spherical harmonics: 



61{n) _ I{h) - (/) ^ 

-=- = = = 2^ aerr,Yern[n). (8) 



im 



Therefore (dimensionless) expansion coefficients aim are 
obtained from the intensity map through the inverse re- 



lation. 



1 



(^)/sky 
1 



sky 



dQf,6I{n)Y;„,{n) 



dnf,I{n)Y,Un), 



(9) 



where /sky = f^sky/47r, and the second equality holds 
except for monopole (i.e., for ^ ^ 0). 
From Eq. ([9|), we have, for ^ > 1, 



•f sky 



(10) 



where C(ni,n2) = (/(ni)/(n2)) /(/} . Thus, we want 
to evaluate (/(ri.i)/(n2)), which now through Eq. ([2]) de- 
pends on 

S^'ixi - x,)5{Li - L,)S^{X2 - Xj)5{L2 - L,)j 
dnsh{xi,Li) dnsh(a^2, ^2) , 



+ 



dL 

dnsh{xi,Li) 
dL 



dL 



-(1+Csh) 



S^{xi - X2)6{Li - L2 



(11) 



Since Xi and X2 represent positions of the subhalo seeds, 
the first term correlates positions and luminosities of two 
distinct subhalos (two-subhalo term). Here (_sh is the 
intrinsic two-point correlation function of the subhalo 
seeds. This two-subhalo term corresponds to a "one-halo 
term" of the halo model which is proportional to 
'T-sh(iCi)nsh(a;2) {N{N — 1)) and gives dominant contri- 
bution to the galaxy power spectrum at scales smaller 
than the virial radius of halos. Although in the halo 
model one often discusses galaxy power spectrum and N 
stands for the total number of galaxies in the host halo, 
exactly the same argument can be applied to the subhalo 
power spectrum, and therefore, TV is regarded as num- 
ber of subhalos instead. The correlation function ^sh is 
related to TV via 1 ^^h = {N{N - 1)) / (iV)^ and nu- 
merical simulations show that 1 -|- ^sh is very close to 1 , 
if the host is massive enough to contain large number of 
galaxies (subhalos). A'' 1 [6l|,[6l|, the case we consider 
here. The second term of Eq. pT|) . on the other hand, 
represents the case of one identical halo where Xi = X2 
and Li — L2 (one-subhalo term). Therefore, one-subhalo 
and two-subhalo terms of (/(rii)/(n2)) are 



Cish(ni,n2) = 



ds2 / dL / d xL 



C2sh(«-l,n2) = 



-==2 f dsi 

167r2(/) J 

dnsh{x,L) 
X — u[sini — X, M) 

X U(S2"2 - X, M), 

(/(ni))(/(n2)) ,^ , , 
=2 (1 + 4shj- 



(12) 
(13) 
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We work on to further simplify Eq. p2)) . First, we 
again take the luminosity function out of a;-integral with 
the index x k, sifii k, S2n2, as its change on subhalo 
scales is moderate. Then, we take Fourier transforms 
of w(rsh,Af): u{k,M). With this, a;-integral gives the 
delta function that collapses one of the fe-integrals, leav- 
ing a common wave number. Finally, we use the small- 
separation approximation (e.g., [63]), where we take s, ry, 
0, n, and instead of si, S2, mi, and n2; these quantities 
are related via s = (si -|-S2)/2, rj = S2 — si, cos 6* = ni ■n2 
and sifii — S2M2 = r]h ~\- s96. This small-angle approxi- 
mation is valid because, for one-subhalo term, we corre- 
late two points in identical subhalos, and their angular 
diameter is very small, thus 6 <^ I. We now have 



Cish(e,n) 



1 



167r2(/) 



o ds drj / dLL 



_ 2 dn shjsn, L) 
~dL 



(27r)3 



\u{k,M)\'. (14) 



Angular power spectrum is given by 



(15) 



where in the first equality, we decomposed the wave num- 
ber k by the components parallel and perpendicular to 
n, i.e.. A; = fc|| + k±, and used d'^k = dki\d'^k±. 

To summarize, combining Eq. (|16p with Eqs. (|14p and 
(fT7| for the one-subhalo term, and with Eq. (fT3| for 
the two-subhalo term, the angular power spectrum of 
gamma-ray background from galactic subhalos is 



Cf — C, 



Ish 



a 



2sh 



Ish 



167r2/sky(/) 



dL 



T2 



dn 



L) 



dL 



u \ -,M 



(18) 



(19) 



2sh 



X Pi{hi ■ ri2), 



{I) 



(20) 



where in Eq. (jl9p . we again performed solid-angle inte- 
gral first, and used angle-averaged luminosity function 
dush/dL in the integrand. We do not try to further sim- 
plify Eq. ([^0)1 , as the two-subhalo term would be more im- 
portant at large angular scales (as shown below), where 
the small-angle approximation is no longer valid. 



Using Eq. pOj) as well as a relation between spherical 
harmonics and Legendre polynomials, i.e., Eq. (46.7) of 
Ref. [i^l, we have 



1 



sky 



sky 



dfln, J dnn2C{hi,h2)Pi{ni • n2) 
dflf, J d02TreC{e,n)Jo{e9) 

d^eC{0,n)e-''-^. (16) 



sky 

In the second equality, we used an approximation 
Ff(cos0) w Jo(id) valid for small 9 and large ^'s, where 
Pg is the Legendre polynomials and Jq is the Bessel func- 
tion of zeroth order. In the last equality, an gle i ntegral 
between f. and was recovered (see, e.g., Ref. |63|). This 
approximation is again particularly good for small angu- 
lar scales (large ^'s), where the one-subhalo term would 
dominate. Now, it is possible to simplify the 0-integral 
of CishS"^^'^; with Eq. (fT4|l . its relevant part is 



d^9e 



dij 



d^k 
(2^ 



\u{k,M)\'e"'- 



2 ik-{r]n+s0 



dH I drj I ^^^|it(fc,M)|2e^'=ii''e'«-(^'^--^) 



(17) 



(27r)3 

dfc|| / d^k±\u{k, M)\^S{kn)S^{skj_ ^ £) 



u \ k = -,M 



III. RESULTS FOR POINT-LIKE SUBHALOS 

In this and subsequent sections, we apply the formu- 
lae for angular power spectrum derived in the previous 
section to several subhalo models. Here, first, we con- 
sider simple models in which we regard all the subha- 
los as gamma-ray point sources, i.e., u{x,M) = 5^{x). 
Its Fourier transform is therefore u{k,M) = 1, indepen- 
dently of wave number and mass. Then the one-subhalo 
term of the angular power spectrum [Eq. (|19p ] is inde- 
pendent of multipole £, and reduces to the Poisson (shot) 
noise. 



A. Models 

Following results of recent numerical simulations, we 
adopt the power-law mass function, dush/dM cx M~", 
and assume that mass distribution is independent of sub- 
halo positions in the host; i.e.. 



dnsh{r,M) 
dM 



nsh{r) 



1 



M 



(21) 



where we assumed M^ain ^ -^^max and typical value for 
a is about 2. For the subhalo number density nsh{r), 
we adopt two different models. One is an "unbiased" 
model where the subhalo distribution follows the NFW 
density profile of the parent halo, nsh(r) oc yOMW,NFw(''^)- 
The other is an "anti-biased" model where the distribu- 
tion is flatter than NFW profile and features a central 
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FIG. 1: Subhalo number densities for unbiased and anti- 



biased distributions, where a = 2, Mn 



10 Mq, and / = 0.5. Densities as a function of galactocen- 
tric radius, nsh{r), are shown as dotted curves, and angle- 
averaged densities as a function of distance from Earth s, 
riih(s) are shown as solid curves. 



FIG. 2: Gamma-ray intensity as a function of angle from the 
galactic center tp, for both unbiased (solid) and anti-biased 
(dashed) models. Subhalo luminosity is related to its mass via 
L = 10^^ (M//M0) s"^ (/3 = 1), and the other parameters are 
Mmin = W'^Mq, Mmax = 1O^°M0, aud a = 2. The isotropic 
intensity measured with EGRET (for i5 = 10 GeV) is shown 
for comparison (dotted). 



core. For this model, we adopt the Einasto profile [63 |. 
nshir) oc pEin(»') with parameters r_2 = 0.81r2oo,MW and 
QfE = 0.68, where r2oo,MW is a radius within which the 
average density is 200 times the critical density (see Ap- 
pendix The latter model takes into account the ef- 
fect of gravitational tidal disruption of subhalos that is 
stronger and reduces the number of subhalos towards the 
central regions of the halo, and indeed is more likely to be 
the case according to numerical simulations [i^ (see also 
Refs. [6l,[63])- Concrete expressions for nsh{r) for each 
case are summarized in Appendix [Al Figure [1] shows the 
subhalo number densities Ushir), for both the unbiased 
and anti-biased distributions (dotted curves). 

Figure[l]also shows the angle-averaged number density 
nsh(s) (solid). Note that unlike Ushir) that is a function 
of galactocentric radius r, n^(s) is a function of dis- 
tance from Earth s. We adopt the same sky region over 
which we take the average as in Ref. [gl] , where the anal- 
ysis of the isotropic diffuse gamma-ray background was 
performed for the Energetic Gamma-Ray Experimental 
Telescope (EGRET); i.e., the galactic center (|6| < 30° 
and |Z| < 40°, where b and / are the galactic latitude 
and longitude, respectively) as well as the galactic plane 
(|6| < 10°) are masked, for which /sky = 0.75. As the 
galactic center is masked, n^(s) is fairly flat within ~10 
kpc, beyond which it reaches regions where the density 
drops more rapidly pmw oc r~^, and it eventually be- 
comes almost the same as nsh{r). 



We assume positive correlation between gamma-ray 
luminosity and mass, L oc Mf^ {13 > 0). Then, 
the luminosity function is obtained by dns\i/dL = 
{dnsh/dM)\dM/dL\, which yields 

dn^^jr^L) / L yi^-")/f^-^ ^ ^^^^ 

The absolute value of the luminosity, e.g., can 
be kept arbitrary, because we are here interested in in- 
tensity fluctuation divided by the mean intensity, where 
the absolute value cancel out.^ Still we comment that 
in order to make the mean intensity (/) as large as 
the observed value around 10 GeV with EGRET [68j , 
Lmin/M^in has to be no smaller than ~10^^/^^Mq^ 
s~^ if /3 — 1 for the unbiased model [H, |5l|; for anti- 
biased model, on the other hand, this value has to be 
^7 times larger. In Fig. [51 we show (/(«)) as a function 
of V, in the case of L/M = lO'^^f-^M'^ s'^ (/3 = 1) 
for both the unbiased (solid) and anti-biased (dashed) 
models. For comparison, we also show the gamma-ray 
background intensity measured with EGRET for 10-GeV 
photons (dotted) [68| . The number of detectable subha- 



There is also an implicit dependence on it through the lower limit 
of s-integral, s»(L), but it is very weak. 
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FIG. 3: Angular power spectrum of the gamma-ray back- 
ground, where subhalos are assumed to be point sources. 
Solid curves are for the one-subhalo terms Cj"'', dashed for the 
two-subhalo terms Cf^*^, and dotted for the host-halo domi- 
nated case. Thick (red) and thin (blue) solid/dashed curves 
are for the unbiased and anti-biased subhalo distributions, 
respectively. Top and bottom solid curves correspond to the 
cases of Mmin ~ W^Mq and 10~^ Mq, respectively, while the 
other parameters are common M^ax = 10^° Mq, a = 2, and 
P = l. 



los and angle-averaged mean intensities for these models 
are (TV^h) = 100, "(7) = 6.4 x 10"^ cm^^ sr"! for the 
unbiased, and (TV^h) = 12, "(7) = 8.7 x 10"^ cm-^ s'^ 
sr^^ for the anti-biased distribution. 



B. Angular power spectrum 

We now move on to the angular power spectrum, start- 
ing with the two-subhalo term [Eq. (|20p ]. which would 
dominate for large angular scales. To evaluate Cf^^, 
we use HEALPix^ package [6^ to generate and ana- 
lyze the gamma-ray map, for which we used parameters 
Nsidc = 1024 and Npi^ = 12iVsMo - 1-2 x lO"^ that corre- 
spond to pixel size of 0.057°, small enough compared with 
angular resolution of Fermi, ~0.1° for a 10-GeV photon. 
The resulting Cg^^ has a highly oscillatory feature as a 
function off, so we average over 0.5£ logarithmic bin for 
a given £. 

In Fig.H we show i{£ + l)Cf' V27r, for both the unbi- 



^ http://healpix.jpl.nasa.gov 



FIG. 4: The Poisson term of the angular power spectrum 
Cp [Eq. (f23)) ] for unbiased (solid) and anti-biased (dashed) 
distributions. Top and bottom sets of curves are for a — 1.8 
and 2, respectively. For the luminosity- mass relation, L = 
102^(M//Mq) s-1 is adopted. 



ased and anti-biased models (dashed curves), assuming 
^sh ^ 1- As expected, it becomes anisotropic for large 
angular scales and small multipole ranges. At smaller an- 
gular scales, where I'^C^^^ cx £^^, the two-subhalo term 
becomes less important compared with the one-subhalo 
(Poisson) term for which i'^Cj^^ oc £^ as we see below. 
Therefore, in general, the two-subhalo term can be safely 
neglected for small angular scales. The difference be- 
tween the unbiased and anti-biased distributions is not 
very large, a factor of two larger for the former. We also 
confirmed that the dependence on a chosen mask is weak. 

There is also contribution to the gamma-ray back- 
ground from a smoothly distributed dark matter com- 
ponent. Its emissivity profile is then proportional to 
p^-y^(r) as annihilation is a two-body process. The angu- 
lar power spectrum from this smooth component can be 
also evaluated with Eq. ([^0]) . but by using line-of-sight 
integral of p^w ^'^^ (^(''^)); rather than of pMW or pEin as 
we see in, e.g., Eq. (O. Figure [3] also shows the angular 
power spectrum from the host-halo component (dotted) . 
The tendency is the same as the two-subhalo terms, with 
amplitudes further larger by a factor of ^3 than the un- 
biased subhalo distribution. Note that in this case, the 
mean intensity (/) is evaluated assuming that only the 
smooth halo gives contribution to the gamma-ray back- 
ground (i.e., no substructure). If we consider both the 
subhalos and smooth halo, the amplitude is reduced de- 
pending on contribution to the mean intensity from each 
component (see Sec. lIVBj) . 
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We now discuss the one-subhalo term Cj^^. If all the 
subhalos are regarded as point sources, u{k,M) = 1, 
Cj^^ becomes independent of £, Cj^^ = Cp, where the 
Poisson-noise (shot or white-noise) term is given by 



Cp = 



167r2/,ky(/)^ 



dL 



ds „ 



dL 



(23) 

Its integrand depends on subhalo number density rtgh, 
whereas in the denominator (/) appears and each (/) 
depends on the subhalo density. Thus, roughly speak- 
ing, Cp is inversely proportional to number of subha- 
los that give significant contribution to the mean inten- 
sity. In Fig. H we also show l{i + 1)0}''^ /2t: for the 
unbiased (thick-red) and anti-biased (thin-blue) distribu- 
tions. By fixing the parameters Mmax = lO'^^M©, a = 2, 
and /3 = 1, we compare results for Afmin = 10^ Af© and 
10~^Mo. In the case of larger Mmin, the mean intensity 
is dominated by bright (massive), relatively rare subha- 
los. On the other hand, in the case of small Mmin, one 
includes fainter subhalos, which increases effective num- 
ber of sources that contribute to the mean intensity and 
thus reduces Cp. This trend is clearly seen in Fig. [3] and 
consistent with the earlier report in Ref. [43|. In Fig. 21 
we plot Cp as a function of Mmin for the unbiased and 
anti-biased distributions and for a — 2 and 1.8. Here 
we rescaled / such that contribution to it from the mass 
range of lO^-lO^^M© is 0.125; by this, for instance, we 
get / = 0.5 (0.15) for a = 2 (1.8) and M,„i„ = IQ-^Mq. 
When the luminosity function is biased towards high- 
luminosity range as in the case of a = 1.8, including 
smaller subhalos has less significant impact on the mean 
intensity. Therefore, the angular power spectrum is fairly 
flat for Mmin < lOM© in this case. 



IV. RESULTS FOR EXTENDED SUBHALO 
MODELS 

A. Subhalo models 

We move on to evaluating the angular power spectrum 
with more realistic models where angular extension of the 
gamma-ray intensity profile for each subhalo is taken into 
account. The models we use are the same as those given 
in Sec. IIII Al except for the source extension u{k,M) as 
well as the luminosity-mass relation that then affects the 
luminosity function. 

As the subhalos are in gravitational potential well of 
the host, they are subject to tidal disruption, and there- 
fore, their outer regions are stripped away. The inner 
regions, on the other hand, are more resilient against 
such an effect. Hence it would be a good approximation 
to assume that the subhalo density profile is given by a 
truncated NFW profile: 

Psh(rshlAf) = ( ^NFw(rshlM) for < r^uu (24) 
^ ^ ' ^ \ for rsh > rcut, ^ ' 



where rcut is a cutoff radius and this is typically much 
smaller than the virial radius. This has been studied 
extensively in Ref. [z3l (see also Refs. [zllTliTi)- 
The gamma-ray luminosity is then given by 



L 



{av) 



/f2 



d^shpL(r.h|M) 

1-1/(1 + Ceut)' 



247r rf [ln(l + Ccut) - Ccut/(1 + Cc 



(25) 



where (av) is thermally averaged annihilation cross sec- 
tion times relative velocity, is WIMP mass, and iV^ is 
number of gamma-ray photons emitted per annihilation. 
In the second equality we define 



K 



(av) 



(av) 



3 X 10-26 cm3 s-i 



100 GcV 



06 



(26) 

where we normalize each parameter with typical values 
often taken in the literature, and Kq = 1.8 x 10~^° 
cm'^ s^^ GeV-2. The value of {av) is closely related 
to the relic density of dark matter as it determines 
the abundance of dark matter particles having survived 
pair annihilation in the early universe. For the super- 
symmetric neutralino, whose mass is around 100 GeV, 
(av) = 3 X 10~26 cni^ g-i jg tj^e canonical value Jj, while 
a wide range of parameter space is still allowed [3]. To 
obtain = 0.6, we used the result of Ref. [11], where 
the gamma-ray spectrum as a result of hadronization and 
decay of 7r° has been fitted with a simple formula. We 
also introduce additional boost factor bsh due to internal 
structure in the subhalo (see below). 

The luminosity also depends on the volume integral of 
the subhalo density squared Pgjj(rsh), which is rewritten 
in terms of mass M, scale radius Ts, and "concentra- 
tion" parameter Ccut = fcnt/rs of subhalos.^ Both rg 
and Ccut are functions of subhalo mass. In order to ob- 
tain values of these quantities for a given mass, we adopt 
scaling relations among various quantities found in the re- 
cent numerical simulations of Ref. (65j . More specifically, 
we adopt the following empirical relations between VJnax, 
r^ax and M: M = 3.37 x 10^MQ(Knax/10 km s-'^ f-^^ 
and (l^nax/ifow)' = 2.9xl04(M/108 Mo)-" !^ where 
VJnax is the maximum rotation velocity of the subhalo, 
?'max is the radius at which the rotation curve hits the 
maximum, and Hq is the Hubble constant. We here 
postulate that for most of the subhalos investigated in 
Ref. [6^ the density profile is well approximated by NFW 
within Tmax (i-C-, '"cut > ''max) and it is indeed the case 



^ We note that concentration parameter is conventionally defined 
as a ratio of virial radius and scale radius. Here, we use the same 
terminology also for a different, albeit similar, quantity Ccut- 



for a sample of subhalos shown in Fig. 22 of Ref. [65j . 

Therefore, we can use the relations between (l^axj ''max) 
and {ps,rs) for the NFW profile, Vg = rniax/2.163 and 
Ps = (4.625/47rG)(V„iax/''s)^, to obtain ps and as a 
function of subhalo mass M. Then finally, to get the 
cutoff radius rcut, we require that the volume integral of 
Eq. (I24p equals to M. We find that the concentration 
parameter Ccut is a decreasing function of M, but still 
larger than 2.163 (i.e., rcut > 2.163rs = r,„ax) even at 
the resolution limit of the simulation, AI — A x IO^Mq, 
confirming that this procedure gives consistent values for 
Ps, rs, and Ccut- These empirical relations do not hold 
at lower mass regions. Therefore, we adopt two different 
approaches, one simply extrapolating the same relations 
down to Earth-mass scale 10~^Mq, and the other assum- 
ing there is no contribution from subhalos less massive 
than IO^Mq. 

Recent numerical simulations also tend to imply pres- 
ence of internal structure within subhalos, i.e., "sub- 
subhalos." If we include these sub-subhalos, both the 
subhalo luminosity and its spatial profile will change. For 
the former, it will give additional boost bsh, for which we 
adopt bsh — 5 (2) for a — 2 (1.9 and 1.8) as it is weakly 
dependent on subhalo masses [l^l- Hence, this effect on 
the angular power spectrum is expected to be rather mi- 
nor. For the latter, on the other hand, the effect would 
be more prominent. If the luminosity is dominated by 
a sub-subhalo component, it would change the smooth 
density-squared intensity profile u{rs\nM) oc Psij(rsh) to 
the clumpy one that is close to proportionality with the 
density, i.e., u(r'sh, M) cc Psh(?'sh)- Although tidal disrup- 
tion is likely to convert the sub-subhalo component into 
the smooth subhalo component without much changing 
luminosity from the latter [tI], we still adopt this model 
regarding it as an extreme possibility. The Fourier trans- 
forms of the emissivity profile for both cases are summa- 
rized in Appendix [Bl 

Considering all of the recent development summarized 
above, we here adopt as our fiducial subhalo model 
the following characteristics: (i) power-law mass func- 
tion dush/dM oc M~^'^ with maximum mass of 
-Mmax = lO^^M©; (ii) anti-biased subhalo distribution 
in the galactic halo, Ush oc pEin; (hi) scale radius rg and 
concentration parameter Ccut for subhalos as a function 
of mass obtained from the empirical relations in Ref. [65j ; 
(iv) subhalo luminosity obtained with canonical values of 
the particle-physics parameters, K — Kq, and the addi- 
tional subhalo boost of fesh = 2; and (v) emissivity pro- 
file of each subhalo dominated by its smooth component, 
u oc Pgjj. We refer to this fiducial set of parameters as 
the model "A." However, we still do not know if the phe- 
nomenological relations visited in (iii) above as well as the 
mass function still hold at the mass scales smaller than 
the current resolution limit. Therefore, we adopt two 
minimum mass scales: Mmin = 10~^Mq for the model 
Al and Mmin = IO^Mq for the model A2. In Sec. HVBl 
we extensively discuss the angular power spectrum for 
these fiducial models, and in Sec. llVCl we study depen- 
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FIG. 5; Luminosity-weighted mass function for subhalos, 
Ldnsh/dln M at s = 0, following the results of the numer- 
ical simulation [6^ (S08). The models based on calibrations 
of field halos, Refs. [7a | 
comparison. 



(BOl) and n§\ (N07) are shown for 



dence of results on chosen models and parameter values. 

Before closing this subsection, in Fig. [51 we show 
the luminosity- weighted mass function, Ldn^/d In M at 
s = obtained for the fiducial models. This quantity 
tells us which mass range contributes to the gamma-ray 
intensity the most [see Eq. ©J. The fiducial subhalo 
model is labeled as S08 according to the reference on 
which this model is based (Ref. [6a]), and this function 
is well fitted by the M'^-^^ scaling (L oc M" "^^). There- 
fore, the smaller subhalos give more important contribu- 
tion to the mean intensity. In the same figure, we also 
show results based on other mass-concentration relations 
given in Refs. [zl (BOl; see also Refs. [IQ,!??,]) and (zl 
(N07). These are, however, for field halos that are not 
fell in potential well of another larger halo, and hence, 
their mass is virial mass and the concentration parame- 
ter is defined as the ratio of virial radius and scale radius 
as conventionally done. These models shown in Fig. \S\ 
thus tell us what the luminosity-weighted mass function 
would be i/ there were no tidal forces acting on subhalos. 
Note also that these concentration models are calibrated 
at even larger scales, such as of galaxies and galaxy clus- 
ters, and the results shown here is based on even more 
violent extrapolation. Nevertheless, it is shown that the 
subhalos are more luminous than the field halos of equal 
mass, and this difference might be as large as two or- 
ders of magnitude. Qualitatively, this difference can be 
explained as follows. Tidal force strips the outer region 
of subhalos away, but the central region is more strongly 
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FIG. 6: Intensity spectrum of the gamma-ray background for 
subhalo models Al and A2, compared with the EGRET data. 
These models are boosted by K/Kq = 13 (Al) and 53 (A2). 



bound. This will reduce the mass of the subhalos signifi- 
cantly but hardly affect the gamma-ray luminosity that is 
proportional to the density squared. Finally, in Fi g.[6 l we 
show energy spectrum of the mean intensity E'^d{I)/dE 
for the models Al and A2, compared with the EGRET 
data [68] . These subhalo models are boosted by a factor 
of K/Kq = 13 (Al) and 53 (A2), with which associated 
anisotropics would be detected (see discussion in the next 
subsection) . 



B. Results for the fiducial model and detectability 
with Fermi 

In Fig. Ula.), we show i{£ + 1)Ci/2tt for the fiducial 
model Al. The two-subhalo term [Eq. (gOl) with ^sh < 1] 
is much smaller than the one-subhalo term [Eq. p9)) ] for 
large multipole ranges. For comparison, we also show the 
Poisson noise [Eq. evaluated for the same model, 

which would be realized if all the subhalos were to be 
gamma-ray point sources. As expected, the power spec- 
trum is more suppressed at smaller angular scales (higher 
multipoles) compared with the noise-like spectrum. This 
means that internal structure of the subhalos should be 
probed with this analysis. 

In fact, we can understand this qualitatively, by an- 
alyzing the integrand of Eq. ([TO]) . In Fig. [51 we show 
contributions to Cj^^ from unit logarithmic mass range 
and from unit logarithmic distance (s) range. The mass 
distributions [Fig. [5)Ja)] peak at high-mass range close to 
Afmax, but are broader for smaller angular scales. This 



FIG. 7: (a) Angular power spectrum for the fiducial subhalo 
model with Mmin = 10"® M© (Al of Table IJ). Contribu- 
tions from the one-subhalo and two-subhalo terms are shown 
as solid and dashed curves, respectively, while dotted curve 
shows Poisson noise that would be obtained if subhalos were 
point sources, (b) Errors for the angular power spectrum of 
the signal SCyCl, for /^h = 0.5 and /b = 0.5. The hori- 
zontal arrow represents the bin width (A^ = 0.5£) for error 
estimates. 



is because at small angular scales, massive subhalos are 
regarded as extended, suppressing the power; note that 
\u{l/s^M)\'^ is a decreasing function of M for fixed l/s. 
Subhalo masses averaged over this distribution and cor- 
responding scale radii are 1.5 x lO^M© and Ts = 1.5 kpc 
{t ^ 10), 1.2 X IO^Mq and r, = 1.4 kpc {i = 100), 
and 6.4 x IO^Mq and = 1.1 kpc {£ ^ 1000). Now, 
Fig. ^h) shows that the contribution from farther sub- 
halos is more important for smaller angular scales, since 
the closer subhalos are more extended. Features at 15 kpc 
correspond to s*(ii„ax), below which contribution from 
massive subhalos are not included as they are identified 
as individual sources. Distances averaged over this distri- 
bution are s = 13 kpc {i = 10), 20 kpc {£ = 100), and 32 
kpc {£ — 1000). Combining these typical distance scales 
with the scale radii, we find that the angular extension 
of the subhalos is typically 6.6° {£ = 10), 3.9° (£ = 100), 
and 1.9° {£ = 1000). For the latter two scales, the sub- 
halo extensions are larger than the angular scales probed 
{6 « 180°/^) and thus typical subhalos are extended, but 
for the case oi £ ~ 10, they are almost point-like sources. 
Therefore, as we see in Fig. [Tja), the one-subhalo term 
starts to deviate from the white noise above £ ^ 10. 

In Figs. IHl^a) andlTOl we show the angular power spec- 
trum, and mass and radius distributions, respectively, 
for the other fiducial model A2 (Mmin = lO^'M©). The 
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FIG. 8: Contributions to from (a) unit logarithmic mass 
range, and (b) unit logarithmic distance range, for the fiducial 
subhalo model Al. The curves are for t — 10, 100, and 1000. 



amplitude of the angular power spectrum for the one- 
subhalo term is much larger than that for the model Al, 
whereas the spectrum shape is almost unchanged. This 
dependence and its interpretation are the same as those 
discussed in Sec. HUB] for simplified subhalo models (see 
Figs. [3] and [4]). The mass and distance distributions for 
Cj^^ are almost the same as the case of the model Al. 

We now discuss the detectability of the angular power 
spectrum. The one-sigma errors of Ci can be estimated 



FIG. 9: The same as Fig.[7]but for the fiducial model A2 with 
Mmin = 10'' Mo, and (b) /sh = 0.1 and /b = 0.9. 



as 



SCf 



{2£ + l)A£Uy 



+ 7772 



(27) 



where is the bin width for that we take = 0.5^, 
and Wt = exp(— o'^£^/2) is the Gaussian window func- 
tion with the beam size Cb — 0.1° of the Fermi (for 
10-GeV photons). The noise power spectrum Cn is 
associated with the finite photon count and it is re- 
lated to number of signal {Ns) and background photons 
(TVbg) through Cn = (aky/^s)(l + ^bg/^s). For 
we assume that half of the EGRET gamma-ray back- 
ground intensity remains unresolved with Fermi, thus 

Ns = /pormiAffteff^^sky « 1.4 X 10^, where /Fermi > 

10 GeV) = 0.5/egret = 0.5 x 10"^ cm'^ s'^ sr-^, 
~ 10''' cm^ is effective area, icff — T'f^fov/47r k, 3x 10*" 
s is effective exposure time, for which we assume 5- 
yr all-sky survey (T = 5 yr and f2tov = 2.4 sr is the 
field of view of Fermi). The background due to detec- 
tor noise is negligible for Fermi, and at high latitudes, 
the galactic foreground due to cosmic-ray propagation 
will be relatively small compared with the isotropic com- 
ponent of the gamma-ray background [7§| (but see also 
Ref. [13|)- Hence we approximate iVbg ^ Ng and obtain 
Cn = 6.7 X 10-5 sr. 

We here consider a multiple-component scenario in 
which the observed total gamma-ray intensity comes 
mainly from several origins. We assume these are dark 
matter subhalos, smooth dark matter component in the 
host, and an astrophysical source such as blazars: 



'host 



(28) 
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FIG. 10: The same as Fig. El but for the model A2. 



and we define eacfi fraction by 

p -^sh p -^host p -^b ( 

Jsh = y J /host = ~y : Jh = T j l^"J 
Jtot Jtot Jtot 

where for simplicity, we represent tiie angle and 
ensemble-averaged intensity (/) by /, and the subscript 
"b" stands for blazars. Then, the total angular power 
spectrum corresponding to the total intensity /tot is given 
by 

Cr = flCt + /LtC^' + /bQ' + (cross terms). (30) 



As galactic dark matter and blazars (or any other extra- 
galactic sources) are independently distributed in space, 
cross-correlation terms, 2fshfhCf^'^ and 2fhostfhC^°^*^''^ 
can be safely neglected. The other cross term 
2/host/shC'|'^'^°''* will be nonzero, but the amplitude of 
^sh.host ^Yionld be at most comparable to C^""^ or Cf 
Therefore, we shall neglect all the cross-correlation terms 
in the following discussions; this is conservative because 
adding this additional component would work favorably 
for dark matter detection. Cf^ is the subhalo angular 
power spectrum and is the same as Eq. (jlSp that we 
closely investigated thus far. The host-halo intensity /host 
is given by 



-^host — 



il-f?K 
Svr 



dsp'^ 



MW 



(s), (31) 



and for the models Al and A2, we have /host//sii ~ 
6.4 X 10~^ and 0.16, respectively. Therefore, the host- 
halo component would be small in the mean intensity, 
and therefore would be further suppressed in the an- 
gular power spectrum. To estimate the blazar power 
spectrum C^, we use the luminosity dependent den- 
sity evolution model for the gamma-ray luminosity func- 
tion [sil, [111, and approximate as nearly a Poisson 
noise, ~ 2 x lO^'^ sr, which is good especially for 
e > 30 l^. We should use the Cf* in the right-hand 
side of Eq. ((27|) to estimate the errors for subhalo signal 
6Cf^. Here we are not interested in the blazars and treat 
them as a reducible component by using lower-energy 
data; see Ref. [i^ for a more detailed discussion. 

In the following analysis, we take fsh as a free param- 
eter instead of others such as fogh and K. We neglect the 
host-halo term as it is always smaller than the subhalo 
term, and so /b = 1 — fsh- In Fig. [TJ^b), we show errors of 
the angular power spectrum, SCf^/f^Y^Cf^, for the model 
Al, when fsh — 0.5 and /b — 0.5. This shows that if 
subhalos give a fractional contribution to the mean back- 
ground intensity, it could also be detected in the angular 
power spectrum, especially at ^ ~ 10. Figure [9l^b) is the 
same but for the model A2 and fsh = 0.1 and /b = 0.9. 
In this case, the detection is more promising as the sub- 
halo anisotropy is much larger than that of blazars, and 
signal-to-noise ratio exceeds 1 for a wide multipole range, 
reaching maximum at ^ ~ 100. Therefore, we define de- 
tection criterion by setting — 1 at either ^ = 10 or 100, 
where ai = f^-^Cf^/SCf^ is the signal-to-noise ratio, and 
for the fiducial models, the value of fsh required to satisfy 
this criterion is fsh = /.t* = 0.24 (Al) and 0.038 (A2). 
(Although it is for only one-sigma detection, we could 
also use other multipole bins to increase significance.) 
For the model Al, in order to achieve fsh — 0.24, we 
need additional boost of /CZ/Cq = /Cdot/^o = 13. For this 
boost factor, the expected number of subhalo detection is 
(iVsh(/fdct)) =0.64. We also find that associated gamma- 
ray flux from t he g alactic center still satisfies constraint 
from EGRET [11], i.e., FGciKdct)/F§S^^'^ = 0.028, 
even if the r^^ cusp of the NEW profile extends to the 
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very center. Thus, the angular power spectrum could 
be a stronger probe than detection as single identified 
sources. For the model A2, the values for these quantities 
associated with the anisotropy detection are K^ct/K^ — 
53, (Ar,h(ifdct)) = 4.0, and FGc(ifdot)/F|g^ET ^ o.l2. 
As this model features smaller mean intensity, we need 
relatively large boost to give a small fraction (~4%) for 
anisotropy detection. Accordingly, the associated num- 
ber of subhalo detection is a few, but it is still not very 
many. Furthermore, the number of subhalo detection po- 
tentially fluctuates to give N^h = 0. Even in this case, the 
anisotropy analysis, therefore, might provide equally sen- 
sitive, but statistically more stable method to probe dark 
matter annihilation in the galactic substructure. The 
spectra of the mean intensity for these models Al and 
A2 with the boost of Kdct/Ko are shown in Fig. [SI 



C. Dependence on models and parameters 

In this subsection, we investigate dependence of results 
on models and parameters for subhalos. In Table [H we 
show models with that we investigate the dependence 
on a, subhalo distribution (unbiased versus anti-biased), 
presence of sub-substructure, as listed in the second to 
seventh columns of the table (we fix Mmax — IQ^'^Mq). 
The first two models are the fiducial models, on which 
we focused in the previous subsection. The next two "B" 
models are for the unbiased subhalo distribution. Models 
Cla and Clb is the same as Al but for different mass 
function. Lastly, the model D2 is the same as A2 but the 
subhalos feature much more extended emissivity profile, 
proportional to the density psh, which might be the case 
if there are a lot of sub-substructure remaining in the 
subhalos. 

In Fig. 1111 we show angular power spectrum (one- 
subhalo term C}^^) of Bl and B2 (unbiased), compared 
with Al and A2 (anti-biased). For all the models in- 
vestigated, the angular power spectrum is larger for the 
anti-biased distribution than the unbiased one. This ten- 
dency is the same as we have seen in the previous section 
for point-like subhalos. Furthermore, deviation from the 
white noise at small angular scales is slightly more sig- 
nificant for the unbiased distribution. This is because 
the relevant subhalos are closer in the case of the un- 
biased model (see Fig. [T|), and therefore, they are more 
extended. 

In Fig. [T^l we compare the models Cla and Clb with 
the fiducial model Al, to study dependence on a. The 
shape of the power spectrum is almost the same among 
these models and the amplitudes follow the same ten- 
dency as seen in Fig. [H Figure [T3] shows a more ex- 
tended model D2 where we assume that the gamma-ray 
luminosity is proportional to the density psh- The sup- 
pression of the power spectrum at small angular scales 
is even more prominent for the more extended emissivity 
profile, which makes it more difficult for this model to be 
detected with anisotropy signals. 
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FIG. 11; Angular power spectrum (one-subhalo term; Cj'^^) of 
gamma-ray background for the models Bl and B2, compared 
witli Al and A2 (Table|l]). The two-subhalo term is also shown 
as dotted curves for both anti-biased and unbiased subhalo 
distributions. 

Table U also shows results of host-to-subhalo ratio for 
the mean intensity, /host//sh- In any of these models, the 
host-halo component is always smaller than the subhalo 
component in the mean intensity. It is, therefore, even 
more suppressed in the angular power spectrum. We also 
show values of f^^^ and the particle-physics parameter 
^dot (in units of Kq) necessary to boost the subhalo sig- 
nals to the level of ct^ = 1 where £ is either 10 or 100. For 
any of these models, we find that if the subhalo contribu- 
tion is as high as ~20%, it should be detected in the angu- 
lar power spectrum, a similar conclusion as obtained for 
the extragalactic dark matter signal [i^ . For the models 
that feature large Cf^ such as A2, B2, and Clb, the sen- 
sitivity of Fermi for fsh could reach down to several to 
ten percent level. We note that to boost the anisotropy 
signal to the detectable level, we only need modest val- 
ues for the particle physics parameter, K^ct/KQ ~ 1~ 
100. The quantities {N^h) and Fqc/F^^^^'^ associated 
with anisotropy detection are summarized in the last two 
columns of Table HI The former ranges from less than 
one to ~10. All the models satisfy a constraint from the 
gamma-ray flux from the galactic center. 



V. DISCUSSION 

Our results confirm what was found in the previous 
papers [13,14^; i.e., the angular power spectrum is dom- 
inated by the noise-like term with a stronger suppression 
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TABLE L Subhalo models considered in calculations of the angular power spectrum. The first column represents (1) model 
identification name. The fiducial models are Al and A2. The second to seventh columns specify each model: (2) value of 
Afmin, (3) a, (4) /, (5) whether subhalo distribution is unbiased or anti-biased {ush oc pnpw or pEin), (6) whether sub-subhalos 
dominate luminosity or not (u oc or p^h), and (7) the boost factor for subhalos 6sh- The eighth column shows the values 
of (8) /host//sh. The rest is the values of (9) f^^\ (10) K^,t/Ko, (11) {N,i,{K^,t)), and (12) FGc{K^ct)/F§§^'''^ , when the 
subhalo contribution is detected in the angular power spectrum (i.e., at — fshCf^/SCf^ — 1) at either £ = 10 or 100. 



Model Mniin 


a 


/ 






) 6sh 


/host 
/sh 


rdct 
Jsh 


Kg 


(Afsh(A'dct)> 


f'GcK-Ffdot) 


Al (fiducial) 10"'' M© 


1.9 


0.2 


PEin 


2 
Psh 


2 


0.0064 


0.24 


13 


0.64 


0.028 


A2 (fiducial) lO^M© 


1.9 


0.16 


PEin 


2 
Psh 


2 


0.16 


0.038 


53 


4.0 


0.12 


Bl IO-'^Mq 


1.9 


0.2 


Pnfw 


Psh 


2 


7.9 X 10"" 


0.21 


1.4 


0.72 


0.010 


B2 10" M© 


1.9 


0.16 


Pnfw 


2 
Psh 


2 


0.021 


0.069 


13 


12 


0.033 


Cla lO-'^M© 


2.0 
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PEin 


2 
Psh 


5 


6.9 X 10"^ 


0.24 


0.37 


0.016 


0.0014 


Clb lO-'^M© 


1.8 


0.15 


pEin 


2 
Psh 


2 


0.068 


0.090 


47 


3.2 


0.11 


D2 10* M© 
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0.16 


PEin 
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2 


0.16 


0.11 
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0.39 
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FIG. 12: Angular power spectrum for the models Cla and 
Clb, compared with Al, to show dependence on a. 



FIG. 13: Angular power spectrum for the model D2 compared 
with A2. 



at higher muhipoles. In particular, behaviors of Ci for 
the models where subhalos luminosity profile is propor- 
tional to its density squared^axe very similar to the power 
spectrum obtained in Refs. p7l.l48j. This is quite natural 
because these authors also assumed similar profiles. Our 
formulation, not only gives physical explanations of these 
results, but also enables to investigate other cases such 
that the luminosity profile is simply proportional to den- 
sity, where the suppression at small scales is even more 
pronounced. 

For the flux sensitivity of Fermi, we adopted -Fgcns = 
2 X 10-1" cm-2 s~i for £; > 10 GeV, but this is only 
for sources whose extension is smaller than the beam 



size (point-like sources). If the sources are extended as 
is the case for massive subhalos, then the flux sensitiv- 
ity is worse, and thus s*(I/) becomes smaller. This will 
decrease expected number of subhalo detection, (Ngh), 
because it considerably reduces effective volume. Thus 
the conclusions given in the previous section might be 
rather conservative, and the anisotropy analysis could be 
even better, compared with the detection of subhalos as 
identifled sources. 

Throughout our calculations, we took into account the 
effect of tidal destruction of subhalos by using the anti- 
biased distribution, since stronger gravitational potential 
around a central region of the host halo is expected to 
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work more efficiently to get rid of more subhalos. We 
also took into account the fact that the subhalos feature 
higher concentration than field halos of the same mass. 
This might be more complicated because the tidal force 
is more effective towards the galactic center, and there- 
fore, the subhalos are more concentrated there. This is 
indeed confirmed by recent numerical simulations, where 
concentration parameters scale as a function of galac- 
tocentric radius as oc r"" "^ [83|. Although this effect 
was not taken into account in our calculations, we could 
do that quite easily by using our analytic formulae with 
slight modification. 

We did not consider dark matter annihilation in the 
extragalactic halos in this study, which should also con- 
tribute to the isotropic gamma-ray background. We note, 
however, that this component might be comparable to 
that from the galactic subhalos j39l [i^. being smaller by 
only a factor of a few. The anisotropy structure for this 
component has been investigated in Refs. [4l|, \^ , and 
found that the angular power spectrum could also be 
reasonably large, £{£ + l)C|'^V27r ~ 0.1 for i > 100 [H. 
Including this component further complicates the analy- 
sis, but should be performed when the actual Fermi data 
became accessible. 

Although we focused on gamma-ray photons with E > 
10 GeV, analysis of angular power spectrum can be per- 
formed at any photon energies. In principle, one could 
use an energy spectrum of the angular fluctuation as an- 
other diagnosis [U, [sH, although we do not discuss it 
further. 



VI. CONCLUSIONS 

To conclude, we investigated the angular power spec- 
trum of the gamma-ray background from dark matter 
annihilation in galactic substructure. Our main findings 
are the foUowings. 

1. In contrast to the earlier works [13, ^^Sit re- 
lied mainly on mock gamma-ray maps generated 
numerically, we derived analytic formulae that en- 
able to compute the angular power spectrum di- 
rectly with assumed subhalo models (Sec. [11]). The 
angular power spectrum consists of two terms: one- 
subhalo and two-subhalo terms [Eq. HH])]. 

2. The two-subhalo term Cf' [Eq. ^] depends on 
smooth (ensemble averaged) distribution function 
of subhalos within the parent halo, Ugh- This is 
largely independent of internal density structure of 
the subhalos. We evaluated this term using the 
HEALPix numerical package [g^ and show that the 
contribution is considerably small at small angular 
scales (e.g.. Fig. |3]). 

3. The one-subhalo term Cj^^ [Eq. depends on 
luminosity profile of each subhalo as well as num- 
ber of subhalos that give significant contribution 



to the mean intensity of the gamma-ray back- 
ground. When the subhalo extension can be ne- 
glected (point-like sources), it simply is the Pois- 
son noise Cp [Eq. ([^5]) ]. In this case, if the subhalo 
mass function is close to dush/dM oc M^^ as im- 
plied by recent numerical simulations, Cp depends 
on the lower mass cutoff Mmin as well as whether 
the subhalos follow the unbiased or anti-biased dis- 
tribution in the parent halo (Fig.|4|). 



4. Taking into account radial extension of the subha- 
los emissivity profile suppresses the power spectrum 
at small angular scales, through u(fc, M) (Fourier 
transform of the emissivity profile) of Eq. . As 
fiducial subhalo models, we assume that the emis- 
sivity profile follows internal density squared, as 
well as the M^^-^ mass spectrum, and the anti- 
biased subhalo distribution within the host. We 
adopt 10~^Mq and 10'* Mq as a minimum mass 
of the subhalos. The angular power spectrum is 
suppressed compared with the white noise at an- 
gular scales smaller than ~10° because of the ex- 
tended intensity profile (Figs. [7] and If Mmin — 
1O~^M0, the angular scales around ~10° is favored 
for detection, where the power spectrum is domi- 
nated by the two-subhalo term. For the anisotropy 
detection, the required fractional contribution to 
the mean intensity is ^20%, for which we need 
additional boost of ~10. The associated number 
of subhalo detection as individual sources is less 
than one. If Mmin — 10'^ Mq, on the other hand, 
the amplitude of the angular power spectrum (one- 
subhalo term) is much larger than the former case, 
and therefore, smaller angular scales ~1° would be 
more promising for detection. The fractional in- 
tensity necessary for anisotropy detection could be 
as small as ~4%, for which the boost is ~50 and 
the associated number of subhalo detection is ~4. 
Therefore, the analysis of the angular power spec- 
trum might be stronger than the detection of sub- 
halos as identified gamma-ray sources, and further- 
more, could provide more statistically stable ap- 
proach to the same problem. 



5. We also investigated dependence on subhalo pa- 
rameters by considering several models (including 
the fiducial ones) summarized in Table[Tl We found 
that the amplitude of the angular power spectrum 
is smaller for the unbiased model compared with 
the anti-biased model (Fig. [TT|l . and for a softer 
mass function (Fig. I12p . The angular power spec- 
trum for models featuring more extended luminos- 
ity profile due to dominance by sub-subhalo com- 
ponent is significantly suppressed at small scales 
(Figs. [121) and is more difficult to detect with Fermi. 
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where "f{a,x) is the lower incomplete gamma function, 
and c_2 = ?'vii /^-2- For numerical values of the param- 
eters, we adopt r_2 = 0.8l7-2oo,MW and as = 0.68 [65l |. 



APPENDIX A: NUMBER DENSITY OF 
SUBHALOS 



Here we obtain number density of subhalos by relating 
it to the mass of the Milky- Way halo, for both the unbi- 
ased and anti-biased models. It is obtained by requiring 



/Afvir,MW 



drAnr 



dMM 



M„ 



dM 



(Al) 

and using Eq. ((2T|) . It is also useful to note the following 
relation: 



vir.MW 



^ ln(l 



,Mw) 



Cvir.MW 



1 + C, 



vir.MW 



(A2) 



where Cvir,MW = ?'vir,Mw/''s,MW is the concentration pa- 
rameter for the Milky- Way halo. 

We thus obtain, for the unbiased case. 



, X fPMwjr) 
fishir) = —rj- X 

'J m i n 



(A3) 



for a ^2. 



For the anti-biased case, we assume that the sub- 
halo distribution follows Einasto profile [G^I, PEin{r) oc 
exp[— (2/aE)(?'/r_2)"'^], where r_2 is a scale radius at 
which the density slope is —2. With a proper normaliza- 
tion, we obtain 



/-^^vir,MW 



_3_ 2cT2 



3/aE-l 




APPENDIX B: FOURIER TRANSFORM OF 
EMISSIVITY PROFILE 



The emissivity profile somehow follows density profile 
of subhalos. If it is smooth, then u{rsh,M) oc ^^^(''sh), 
whereas if they include a number of sub-subhalos, then it 
scales more like oc Psh(fsh)- Here we compute the Fourier 
transform of this quantity, u{k,M). When the density 
profile of the subhalos are well described by the NFW 
profile up to a cutoff radius Tcut , then the both cases have 
analytic expressions for u(k, M) given concentration and 
scale radius. These forms are somewhat complicated, and 
thus we instead give fitting formulae that give excellent 
approximation for a wide range of reasonable values of 
Ccut- For the both cases, the fitting form is 



u(fc, M) 



[l + ai(fcr,)i/b + a2(fcr«)2/*']' 



(Bl) 



For u (X (no sub-subhalos), ai = 0.13, a2 = 0, and 
6 = 0.7, and these values are largely independent of Ccut. 
On the other hand, for u oc p (with sub-subhalos), we 



(A4) have ai = 100c~„°-^\ 02 = l7Qc^^^^, and b = 0.16c! 
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